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摘要 : 受 硬件 限制 ,4D-CBCT 成 像 中 单个 呼吸 相位 对 应 投影 数目 过 少 , 而 相应 的 常规 解析 算法 重建 图 像 中 充斥 着 大 量 噪声 和 伪 
影 。 为 解决 此 问题 ,鉴于 当前 呼吸 相位 图 像 可 以 通过 其 它 呼吸 相位 图 像 运动 补偿 获取 ,本 文中 我 们 提出 一 种 利用 多 相位 投影 


图 像 间 的 关联 信息 并 未 引入 新 的 伪 影 和 运动 模糊 。 
关键 词 :四维 锥 形 束 CT; 运 动 补偿 ;压缩 感知 


据 重建 4D-CBCT 的 策略 。 本 文中 我 们 构建 了 包含 基于 多 相位 投影 数据 的 保 真 项 和 基于 压缩 感知 理论 的 全 变 分 正则 化 的 代价 
函数 。 对 于 投影 数据 保 真 项 的 设计 ,不 再 局 限于 当前 相位 投影 数据 ,而 是 利用 多 个 相位 投影 数据 通过 变形 的 成 像 模 型 联合 构 
建 。 对 于 复杂 代价 函数 的 优化 ,我 们 利用 GPBL(Gradient-Projection-Barzilai-Linesearch, GPBL) 算 法 来 实现 。 物 理 体 模 及 临床 
数据 实验 结果 表明 ,相对 于 解析 算法 及 代数 欠 代 全 变 分 约束 算法 ,新 方法 在 噪声 和 伪 影 的 抑制 方面 有 上 佳 表 现 ,引入 不 同 相 位 
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Abstract: Restriction by hardware caused the very low projection number at a single phase for 4-dimensional cone beam 
(4D-CBCT) CT imaging, and reconstruction using conventional reconstruction algorithms is thus constrained by serious streak 


artifacts and noises. To address this problem, we propose an approach to reconstructing 4D-CBCT images with multi-phase 
projections based on the assumption that the image at one phase can be viewed as the motion-compensated image at another 
phase. Specifically, we formulated a cost function using multi-phase projections to construct the fidelity term and the TV 


regularization method. For fidelity term construction, the projection data of the current phase and those at other phases were 
jointly used by reformulating the imaging model. The Gradient-Projection-Barzilai-Line search (GPBL) method was used to 
optimize the complex cost function. Physical phantom and patient data results showed that the proposed approach could 
effectively reduce the noise and artifacts, and the introduction of additional temporal correlation did not introduce new 


artifacts or motion blur. 
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流行 病 学 资料 显示 " ,肺癌 是 我 国 第 一 高 发 的 恶性 
肿瘤 ,肺癌 死亡 率 已 占 城 市 中 癌 亡 的 首位 ,发 病 率 也 比 
10 年 前 翻 了 一 倍 。 肺 瘤 发 病 时 病 期 相对 较 晚 ,放射 治疗 
是 主要 治疗 手段 之 一 ”。 

近年 来 , 随 着 放射 生物 学 .放射 物理 学 的 极速 发 展 


到 直线 加 速 吕 上 ,治疗 中 可 在 治疗 位 置 进行 透视 . 摄 片 、 
容积 成 像 ,继而 在 线 调整 治疗 图 像 与 计划 图 像 的 位 置 ， 
减少 摆 位 误差 ,以 期 实现 精确 放疗 的 目标 ”。 当 前 ,临床 
肿瘤 放射 治疗 多 采用 3D-CBCT 引 导 的 三 维 适 形 放射 治 
疗 技术 。 然 而 ,对 于 胸 腹 部 位 置 的 肿瘤 ,为 实现 精确 放 


图 像 引 导 放 疗 (IGRT) 技 术 的 产生 为 肺 部 肿瘤 精确 放疗 
提供 了 可 能 外 。 目 前 ,IGRT 类 型 较 多 ,以 锥 形 束 CT 
(CBCT) 引 导 的 放疗 系统 为 最 具 代表 性 “, 它 直接 整合 
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疗 , 仅 消除 摆 位 误差 是 不 够 的 ,呼吸 运动 心脏 与 大 血管 
搏动 以 及 食管 蠕动 均 会 引起 肺 部 肿瘤 的 移动 ,其 中 呼吸 
运动 的 影响 最 为 主要 “。 呼 吸 运动 会 导致 重建 CBCT 
图 像 中 出 现 严 重 的 伪 影 及 模糊 ,为 纠正 呼吸 运动 的 影 
啊 , 实 现 自 由 呼吸 状态 下 的 肿瘤 的 实时 监测 ,近年 来 ， 
4D-CBCT 的 概念 被 提出 "，。 然 而 ,基于 现 有 机 载 
CBCT 成 像 系 统 , 因 其 转速 有 限 , 对 于 特定 呼吸 时 相 获 
取 投 影 数目 过 少 , 常 规 解析 重建 算法 尚 无 法 重建 出 高 清 
晰 无 伪 影 的 图 像 来 完成 肺 部 肿瘤 的 准确 定位 及 在 治疗 
时 的 三 维 运动 轨迹 跟踪 (图 1)。 
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图 1 利用 CBCT 旋 转 一 周全 部 投影 数据 重建 所 得 3D-CBCT( 多 
个 呼吸 相位 平均 ) 图 像 ,肿瘤 边界 不 清 ,存在 运动 模糊 和 伪 影 
Fig.l 3D-CBCT image reconstructed from the full projection 


with blurs and artifacts on the edge. 


近年 来 ,为 实现 4D-CBCT 的 优质 重建 ,在 软件 法 
方面 ,主流 的 研究 方法 可 以 归 为 两 大 类 :基于 运动 补 
偿 的 重建 算法 ”和 基于 呼吸 相关 (Respiratory 
Correlated, RC) 的 重建 算法 "中 。 基 于 运动 补偿 的 重 
建 算法 将 4D-CBCT 看 作 一 个 整体 ,利用 全 部 投影 和 
各 呼吸 相位 图 像 间 的 运动 信息 进行 重建 。 其 中 ,由 
投影 数据 或 者 图 像 本 身 估 计 所 得 的 呼吸 运动 矢量 场 
(MVFs) 直 接应 用 于 重建 图 像 ,而 现 有 的 算法 对 其 估 
计 一 般 不 准确 , 易 导致 运动 模糊 校正 不 彻底 或 者 伪 
结构 的 存在 。 另 一 类 为 基于 RC 重建 算法 ,该 类 算法 
利用 对 应 呼吸 相位 的 投影 数据 进行 该 相位 CBCT 图 
像 重建 。 该 类 算法 属于 不 完全 数据 或 者 稀 琉 角度 CT 
重建 范畴 。 针 对 此 类 型 的 问题 ,近年 来 各 种 算法 相 
继 出 现 ,例如 , 随 着 CS 理论 的 发 展 ,Sidky 等 "提出 一 
种 ASD-POCS 算法 用 于 4D-CBCT 重建, 取得 了 不 错 
的 效果 ,条 形 伪 影 得 到 有 效 抑制 ,但 与 此 同时 图 像 出 
现 过 平滑 现象 或 阶梯 效应 。Chen 等 .利用 全 部 投影 
数据 重建 图 像 作为 先 验 图 像 , 提 出 了 PICCS 算 法 ,该 
算法 有 效 的 融入 先 验 图 像 信 息 ,重建 所 得 单 相 位 图 
像 伪 影 得 到 有 效 抑 制 ,但 先 验 图 像 中 的 模糊 伪 影 和 
某 些 虚假 解剖 结构 将 有 可 能 传播 至 待 重建 图 像 , 影 
响 最 终 重建 图 像 质量 。 基 于 呼吸 相关 的 4D-CBCT 重 
建 算法 由 于 投影 数据 严重 缺失 的 本 质 , 导 致 其 重建 
图 像 质 量 不 高 , 且 现 有 的 研究 大 多 局 限于 相对 简单 
的 体 模 研 究 , 对 于 肺 部 复杂 解剖 结构 ,此 类 方法 重建 
图 像 中 细微 结构 有 可 能 会 丢失 或 产生 畸变 。 

本 文中 ,我 们 拟 将 整个 呼吸 周期 的 4D-CBCT 看 
作 整 体 ,充分 考虑 各 呼吸 相位 图 像 之 间 的 相关 特性 ， 
构建 用 于 迭代 重建 的 数据 保 真 模型 ,并 基于 压缩 感 
知 理论 ,引入 全 变 分 正则 化 先 验 ,形成 基于 运动 补偿 
的 压缩 感知 重建 算法 。 文 中 ,我 们 将 利用 物理 体 模 
数据 及 临床 数据 对 算法 进行 性 能 测试 。 
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1 方法 
1.1 基于 运动 补偿 的 压缩 感知 4D-CBCT 重 建 模型 构建 

假设 整个 呼吸 周期 可 分 为 P 个 呼吸 相位 ,临床 上 
该 值 一 般 设 定 为 8-10.4D-CBCT 可 以 定义 为 
f={fi;;…; 记 ee RW ,其 中 入 表示 一 个 呼吸 相位 下 
3D-CBCT 体 数 据 体 素 个 数 . 对 于 单个 呼吸 相位 ,CBCT 
的 成 像 模 型 可 以 用 下 式 表示 : 

ge=Hf,i=1,2,...,P (1) 

式 (1) 中 ,表示 对 应 于 呼吸 相位 ;的 投影 数据 ,所 该 
呼吸 相位 下 对 应 的 系统 矩阵 , /为 待 估计 图 像 。CBCT 
图 像 的 重建 即 由 已 经 获得 的 投影 数据 g; 及 相应 的 系统 
和 矩阵 五, 对 式 (1) 进 行 反 问题 求解 。4D-CBCT 中 个 呼吸 
时 相 下 投影 数据 不 完整 的 特点 ,使 得 式 (1) 的 求解 成 为 
一 个 病态 着 问题 。 

本 文中 ,引入 不 同 呼吸 相位 图 像 间 的 相关 特性 ,我 
们 拟 将 原 CBCT 成 像 模 型 写 为 : 

g,=HM, fo, i=1,2,...P;k=1,2,...P (2) 
其 中 ,假设 f 为 对 应 呼吸 相位 Kk 的 CBCT 图 像 ,利用 相位 
和 j 图 像 间 的 变形 场 ,可 以 得 到 M, ,fi xf;,M 为 图 像 
到 的 运动 补偿 算 子 。 因 此 ,通过 式 (2) 我 们 可 将 不 同 呼 
吸 相 位 的 图 像 进行 关联 ,在 估计 相位 Kk 对 应 CBCT 图 像 
时 ,引入 其 它 相位 对 应 的 投影 数据 以 增加 代价 函数 保 真 
项 的 可 靠 性 ,同时 ,为 保持 图 像 的 边缘 特性 及 去 除 模 
糊 效应 ,我 们 引入 全 变 分 正则 化 的 方法 ,构建 如 下 代价 
函数 : 

minjmize 2 a HM -gi + 让 (3) 
式 (3) 中 ,a 为 权重 因子 ,表征 相位 k 与 相位 图像 的 相 
关 程 度 , 由 于 呼吸 运动 为 周期 性 的 ,因此 ,本 文中 ,与 
呼吸 相位 k 相 临近 的 (如 吸 气 中 期 与 吸 气 末端 ) ,给 予 
较 大 权重 ,与 呼吸 相位 k 相 背离 的 (如 呼 气 末 端 与 吸 气 
末端 ) ,给 予 较 小 的 权重 , 当 i=k 时 aj=1 ,参数 mm 表示 与 
相位 k 相 近邻 的 2m 个 呼吸 相位 。4 为 正则 化 参数 ,用 
来 调节 数据 保 真 项 与 正则 化 项 间 的 权重 比例 。 本 文 
中 ,采用 各 向 同性 全 变 分 来 构建 正则 化 ,基于 TV/L1 
的 正则 化 可 以 表示 为 : 

A  g 
其 中 D, , D,, D, 为 沿 着 横向 纵向 之 轴 方 向 的 前 向 差 
分 投影 算 子 ,定义 为 : 

Df,= vec( fi(x+ 1,y,7) -f(x,y, z)), 

D,fi=vec(f.(x,y+1,7) -f(x,y,7)), (5) 

Df =vec( fi(x,y,z+ 1) —fi(x,y, z)) 
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1.2 基于 光 流 法 的 运动 变形 场 估计 

式 (3) 中 ,不同 相位 图 像 间 运动 变形 场 的 准确 估计 
对 重建 图 像 的 质量 起 着 至 关 重 要 的 作用 。 本 文中 ,我 们 
采用 光 流 法 对 呼吸 运动 进行 估计 。 三 维 光 流 法 
(3D-OFM) 是 一 种 重要 的 分 析 运 动 图像 的 可 变形 图 像 
配 准 方法 , 它 的 计算 是 基于 图 像 像 素 进 行 , 在 精确 性 和 
和 鲁 棒 性 方面 有 很 大 的 优势 。 本 文中 所 用 光 流 法 为 霍 恩 - 
舒 克 (Horn &Schunck) 提 出 的 算法 “。 对 于 一 个 物理 点 
的 位 置 (x, y, zz) 在 一 个 参考 相位 x 和 它 的 位 置 在 任何 其 
他 相位 的 的 对 应 关系 可 用 一 个 变形 模型 定义 : 

f(xtu,yt+v,z+s)=f(x,y, 7) (6) 

(u,v,5) 为 图 像 f 与 f 之 间 的 光 流 场 ,将 待 配 准 的 图 
像 按 照 光 流 进行 变换 即 可 得 到 配 准 图 像 。 转 化 为 算 子 
形式 , 即 为 M ,f=f。 理 论 上 任何 相位 的 CBCT 体 数 
据 f 均 可 由 及 对 应 光 流 场 变 形 配 准 获取 , 且 该 过 程 具 
有 可 逆 性 。 

本 文 提出 的 算法 是 基于 重建 图 像 本 身 进 行 运动 估 
计 / 补 偿 的 ,因此 ,图 像 重 建 和 运动 变形 估计 是 协同 进行 
的 ,在 此 我 们 需要 指出 ,最 初 重建 图 像 中 含有 较 多 的 噪 
声 及 伪 影 ,初始 运动 估计 信息 的 可 能 不 准确 ,但 随 着 迭 
代 次 数 的 增加 ,重建 图 像 质量 的 逐步 提升 ,该 效应 会 逐 
步 减 弱 。 
1.3 算法 优化 求解 

对 于 代价 函数 (3) 式 的 求解 ,由 于 其 复杂 性 及 维度 
比较 高 ,因此 本 文中 我 们 利用 (Gradient-Projection- 
Barzilai-Linesearch, GPBL) 算 法 "进行 迭代 优化 。 对 
于 特定 的 相位 k, 其 优化 具体 过 程 如 下 : 

对 于 k=1,2,.….P, 执 行 

Step (1) : EA 全 > aMi HI (HM, Re -8 


i=k=—m 


tav||DA|, 
Step (2) :对 于 所 有 的 像素 点 j: 如果 满足 
f* "(D>=0,g" (<0, 执 行 pr=gi"" 


否则 p." (j=0 
T 
(n-D) (n—-2) (n) (n-1) 
(4 -fi ) 人 一 有 ) 


2 


Step (3):1" < 


(n—2) 


(2 一 1 
i Sf 
(OO 

8 


Step (Df ef "1 

对 于 所 有 像素 点 广 如 果 满 足 产 "CD<0, 执 行 
A =0 

Step (5): 利 用 f" 及 3D-OFM 算法 重新 估计 M, ,; 反 
复 执行 此 过 程 直至 | ”了 ”省 <e, 或 者 达到 设 定 氨 代 
次 数 。 实 验 过 程 中 ,参数 m 的 取 值 设 定 为 4,4 取 值 范围 


为 1x10'~1x10, 根 据 不 同 实验 而 定 。 


2 结果 
2.1 实验 数据 获取 

本 文中 ,我 们 分 别 利用 物理 体 模 及 真实 临床 数据 进 
行 算法 测试 。 物 理 体 模 设 计 为 静止 的 肺 部 仿真 人 体 体 
模 加 直径 约 为 3 cm 的 “运动 组 织 ” 等 效 肿瘤 。 该 肿瘤 放 
置 在 左 侧 胸腔 ,与 机 械 4D 运动 控制 装置 相连 接 。 肿 瘤 
的 运动 设 定 为 幅 值 为 13.4 m, 沿 轴 向 做 正弦 运动 ,周期 
为 4s。 数 据 的 采集 利用 瓦 里 安放 疗 设备 获得 ,数据 采 
集 参 数 为 : 管 电压 120 kVp, 旋 转角 度 为 360° ,扫描 时 间 
大 约 为 1 min, 共 获取 投影 数据 约 为 620 个 ,原始 投影 妆 
据 大 小 为 1024x768, 像 素 大 小 为 0.388 mmx0.388 mm。 
实验 过 程 中 ,我 们 将 投影 数据 降 维 为 512x384, 相 应 像 
素 大 小 为 0.776 mmx0.776 mm。 利用 标记 点 记录 位 
置 ,我 们 将 所 有 投影 数据 划分 为 8 个 相位 。 

病人 数据 我 们 采用 http://www.openrtk.orgA 网 站 
共享 的 投影 数据 ,扫描 角度 范围 为 360, 总 的 投影 个 数 为 
644, 依 据 网 站 提供 的 投影 数据 相位 信息 ,我 们 将 投影 划 
分 为 10 个 相位 。 

2.2 实验 结果 分 析 

2.2.1 物理 体 模 重建 结果 分 析 图 2 给 出 了 某 一 呼吸 相 
位 的 物理 体 模 的 重建 结果 ,其 中 ,从 FDK 算法 重建 的 图 
像 我 们 可 以 看 到 ,由 于 投影 数据 的 缺失 ,重建 图 像 中 存 
在 大 量 的 伪 影 和 噪声 ,很 多 解剖 结构 信息 无 法 辨识 (图 
2A)。ART-TV 算 法 重建 结果 ,可 以 看 到 条 形 伪 影 得 到 
有 效 的 抑制 ,但 由 于 投影 数据 本 身 的 数据 不 足 ,导致 重 
建 效 果 有 限 ( 图 2B)。 而 本 文 提 出 算法 的 重建 结果 ,可 
以 明显 的 看 出 ,相对 左 侧 图 像 , 右 侧 图 像 噪声 和 伪 影 较 
少 ,结构 清晰 ,图 像 质 量 得 到 很 大 提升 (图 2C)。 进 一 步 
的 , 表 1 给 出 了 以 上 FDK 及 本 文 提 出 重建 算法 在 不 同 
呼吸 相位 下 肿瘤 均 质 区 域 (图 3) 标 准 差 计算 值 ,其 中 所 
选取 得 区 域 为 图 3 所 示 的 区 域 ,大 小 为 20x16 个 像素 。 
从 计算 结果 可 以 看 出 ,本 文 提 出 的 算法 重建 图 像 的 标准 
差 远 低 于 FDK 算 法 重建 图 像 , 进 一 步 说 明 本 文 算法 能 
够 有 效 的 抑制 了 图 像 噪声 。 

图 4 给 出 了 对 应 于 5 个 呼吸 相位 的 本 文 算法 重建 
结果 的 矢 状 面 图 像 (图 4A) 以 及 相应 的 肿瘤 运动 轨迹 图 
(4B)。 结 果 显 示 , 本 文 算法 能 有 很 好 的 重建 出 运动 肿瘤 
而 不 引入 运动 模糊 或 伪 影 ,同时 ,对 比 于 原始 实验 设计 
的 肿瘤 运动 模型 ,本 文 算 法 能 够 准确 的 刻画 肿瘤 运动 的 
轨迹 。 

2.2.2 病人 数据 重建 结果 分 析 图 5 和 图 6 分 别 给 出 了 
病人 数据 5 个 呼吸 相位 的 本 文 算法 重建 结果 的 冠状 面 
和 矢 状 面 图 像 。 图 5 和 图 6 中 ,第 1 行为 FDK 算 法 重建 
图 像 ,第 2 行为 本 文 算法 重建 结果 。 可 以 明显 的 看 出 ， 
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图 2 某 一 呼吸 相位 物理 体 模 重建 结果 


Fig.2 Physical phantom results of a selected phase. A: FDK reconstruction; B: ART-TV reconstruction; C: 


Reconstructions using the proposed algorithm. 


图 3 用 于 计算 噪声 水 平 的 均 质 区 域 ,红色 方 框 表 示 选 定 区 域 
Fig.3 A homogenous region (boxed) used to evaluate the 
noise. 


FDK 重 建 图 像 中 伪 影 和 噪声 掩盖 了 许多 图 像 的 细节 信 
息 , 同 时 使 图 像 的 对 比 度 明显 下 降 ,肿瘤 边界 特征 不 明 
显 。 而 本 文 算法 重建 图 像 , 肺 部 细节 结构 信息 明显 ,图 
像 整体 噪声 水 平 较 低 ,对比 度 较 高 ,肿瘤 边界 相对 明 
显 。 这 对 于 后 续 的 图 像 分 析 ( 如 图 像 分 割 ) 及 肿瘤 运动 
轨迹 提取 都 是 极为 有 利 的 。 


四 
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常规 近代 重建 算法 ”数据 保 真 项 的 设计 一 般 只 
涉及 当前 呼吸 相位 投影 数据 而 忽略 了 不 同 相位 图 像 
之 间 的 相关 性 。 本 文中 ,我 们 充分 利用 4D-CBCT 中 
不 同 呼吸 相位 图 像 间 的 相关 特性 ,通过 引入 运动 补偿 
算 子 ,将 不 同 相 位 的 图 像 关联 在 一 起 ,进而 构建 当前 
投影 数据 与 其 他 相位 图 像 的 关系 式 ,并 将 改进 后 的 成 


图 4 体 模 本 文 算法 重建 结果 的 矢 状 面 图 像 (4) 以 及 相应 的 肿瘤 运动 轨 


Fig.4 Sagittal reconstructed images of the physical phantom (4) and the 
corresponding tumor trajectory curves (B). 


像 模 型 融入 代价 函数 保 丰 项 的 设计 中 ,构建 基于 多 相 
位 投影 数据 的 压缩 感知 4D-CBCT 重 建 模型 。 对 模型 
的 求解 优化 采用 ,GPBL 算 法 实现 。 影 响 算法 性 能 比 
较 关键 的 一 步 为 迭代 重建 过 程 中 的 运动 估计 ,运动 估 
计 不 足 会 造成 重建 图 像 的 运动 模糊 的 存在 ,运动 估计 
过 度 则 会 造成 重建 图 像 出 现 伪 影 伪 像 。 同 时 ,重建 过 
程 中 各 参数 的 选择 也 会 较 大 的 影响 重建 结果 , 参 交 
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图 5 对 应 于 5 个 呼吸 相位 的 本 文 算法 重建 结果 的 冠状 面 图 像 


Fig.5 Coronal reconstructed images of 5 selected phases. The first column shows the FDK reconstructions and the 


second column shows the reconstructions with the proposed algorithm. 


图 6 对 应 于 5 个 呼吸 相位 的 本 文 算法 重建 结果 的 矢 状 面 图 像 
Fig.6 Reconstructed sagittal images of 5 selected phases. The first column shows the FDK reconstructions and the 
second shows the reconstructions with the proposed algorithm. 


表 1 不 同 重建 算法 肿瘤 均 质 区 域 的 标准 差 对 比 


Tab.1 Standard deviations of the homogeneous regions for different reconstruction methods 


Phasestd Pl1 P2 Ps P4 Rs P6 Pp7 P8 
FDK 0.0023 0.0025 0.0023 0.0028 0.0026 0.0025 0.0031 0.002 
Proposed Algorithm 6.62E-4 9.04E-4 9.17E-4 1.00E-3 8.98E-4 8.77E-4 8.62E-4 6.17E-4 


选择 问题 是 一 个 公开 的 难 优化 的 问题 ,在 此 我 们 不 
对 重建 参数 做 过 多 的 讨论 。 本 文 算法 中 ,我 们 引入 
全 变 分 正则 化 的 策略 ,一 则 为 了 抑制 噪声 及 伪 影 ,二 
则 为 了 防止 运动 模糊 产生 ,保持 边缘 特征 。 从 实验 
结果 来 看 ,本文 提 出 的 算法 能 够 有 效 的 抑制 重建 图 
像 中 的 噪声 和 伪 影 ,提高 重建 图 像 的 质量 ,同时 准确 
的 刻画 肿瘤 运动 轨迹 。 
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